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On the basis of the Fermi liquid theory, we investigate the many-body effect on the Drude 
weight. In a lattice system, the Drude weight D is modified by electron-electron interaction 
due to Umklapp processes, while it is not renormalized in a Galilean invariant system. This is 
explained by showing that the effective mass m! for D oc n/m' is defined through the current, 
not velocity, of quasiparticle. It is shown that the inequality D > is required for the stability 
' against the uniform shift of the Fermi surface. The result of perturbation theory applied for 

m 
m 

t^. . around half filling n ~ 1 by Umklapp processes. 

o 
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§1. Introduction 

o 

In comparing theory and experiment of highly correlated electron systems, we must distinguish 
genuine many-body effects from one-body effects. This is trivial at least in a theoretical investiga- 
tion based on a simple model. For example, in the Fermi liquid theory, the thermal mass measured 
in a specific heat experiment reflects the mass of quasiparticles, which can be heavily enhanced by 
interaction effects. On the other side, the effective mass m! defined in the Drude weight, which 
characterizes the metallic property of the system, does not generally coincide with the quasiparticle 
mass m*. In fact it is shown that to' is not renormalized, or it is independent of electron-electron 
interactions in Galilean invariant systems. In a lattice system, Umklapp processes make the Drude 
weight depend on the electron-electron interaction and thus the Drude weight can be used as a 
direct probe of the metal- insulator transition.BB This problem has been studied particularly in 
one dimensional systems ,i'f§i'i) where it is shown that the Drude weight D vanishes as the carrier 
density n approaches half filling, i.e., in the presence of the Coulomb interaction U > 0, we have 
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D oc |1 — n\ = 5 for the low density limit 5 — > of doped hole. Though the metal-insulator transition 
in this form might not be realized in general, we are interested in how this qualitative behavior 
of 'doped insulator' is reconciled with normal Fermi liquids in two or three dimensional systems, 
and we investigate and show how it is put in the context of the general theory of Fermi liquids 
when Umklapp processes are effective. In this respect, the effect of Umklapp processes in Fermi 
liquids was discussed for the cyclotron resonance frequency lo c by Kanki and Yamada,B where it 
was shown that co c is not renormalized in an isotropic system. In a similar manner, we discuss 
that D is not renormalized in an isotropic system and then see in what manner D is modified in a 
lattice system. 

In §|2] we derive the expression for D on the basis of the Landau Fermi liquid theory and show 
that it is calculated as the second derivative of the energy E p of the state \p) , which is defined from 
the ground state by shifting the fermion distribution function n2.- The results obtained in §|| are 
derived microscopically in §|3], where the Fermi liquid result of D is obtained using 



To supplement the general discussion of §§^ and ||, D is estimated explicitly in §||, where the 
effective mass for the Drude weight is calculated up to the second order in U of the Hubbard model 
in a square lattice. We observe that D as a function of n is modified qualitatively around half filling 
n ~ 1. The last section contains discussions. In this article, we investigate a system described by a 
single-band model of fermions at absolute zero. We are mainly interested in the collisionless region 
where the effect of a quasiparticle lifetime can be neglected. In Appendix we outline the derivation 
of one of the main results on the basis of the finite temperature formalism. A general theory at 
finite temperature will be presented in a subsequent paper .B 

§2. Fermi Liquid Theory I 

First we derive the uniform conductivity on the basis of the Landau theory of Fermi liquids.!^ 
In the presence of an applied electric field E, the deviation 5n k from the ground-state distribution 
function 





5n k = n k -n k 



satisfies the Boltzmann equation, 



(q-v k - Lo)5n k + q • v k 5(fi - e k ) ^ f(k, k')5n k > 

k' 

+ieE • v k 5(n - e k ) = 0, 



(2.1) 



where e k is energy of the quasiparticle k and 



de k 
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For q = we obtain 

5n k = *E-v k 6Ui-e k )^ 

u> 



The total current induced by the field is 



J = J2$n k j k . (2.3) 

k 



The conductivity tensor is defined by 



<V(") = r^p ( 2 - 4 ) 



where is the total volume of the system. Hence we find 

, . ie 2 / n \ , 

a A") = — ( — ) > ( 2 - 5 ) 

where 



The Drude weight D, eq. (1^1), is given as the coefficient of the delta function 6 (u) in Reu^^u;); 



n 



m 



Rea^(u) = ne z — (2.7) 



/Ltl/ 



Therefore, we obtain 



D = vre 2 



The effective mass m! , defined by eq. ( |2.6[) , plays an important role in our theory. 
To see that electron-electron interaction does not modify the Drude weight in Galilean invariant 



systems, eq. (|2.q) is transformed as follows, 
Here we used 

5(n-s k )=-^, (2.10) 

i.e., the quasiparticle distribution function n k is represented by the step function 9([i — e k ) at zero 
temperature. As a direct consequence of Galilean invariance, the relation j ku = k u /m, with the 
bare electron mass m, is not modified by electron-electron interactions. Therefore in this case we 
obtain 

—) =-<W (2-11) 
where 

n ^E n °> ( 2 - 12 ) 
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and the Drude weight D = ire^n/m is not affected by the interaction. If we assume 
dj]~u/dk v =const. for simplicity, from eq. (2.9) we find 

-i) (±) = (2.13) 

where m* represents the effective mass of the quasiparticle k. Hence it is concluded that the mass 
ml in the Drude weight is given through the current of quasiparticle, and that m', m* and the 
bare (crystalline) mass m generally take different values. We remark that the effective mass m* 
thus defined by 1/m* oc (v) does not generally coincide with the thermal mass defined in the total 
density of states m£ h oc (v" 1 }, where the average is taken over the Fermi surface. The difference 
between m* and m* h would be significant only in exceptional cases, as in the vicinity of the van 
Hove singularity where the former is less singular than the latter. 

Now let us show that the right-hand side of eq. Q2,6j ) is derived in terms of the state \p), which 
represents the state obtained by displacing the ground state configuration by p in the momentum 
space, i.e., \p) is derived from the ground state |0) by formally replacing the distribution function 
with n\_ p . By construction, the resulting state \p) carries a finite current JpS Denoting the 
energy of the state \p) as E p , we would like to derive J PfJi = dEp/dp^ and express (n/m 1 )^ as the 
second derivative of E p with respect to p. 

According to the Landau Fermi liquid theory, the current carried by the quasiparticle k is given 
byi 

3k = v k + f( k > k')v kl 8(n - e k t). (2.14) 
k' 

The spin dependence of the interaction function f(k, k') was omitted for simplicity. Using 5n k = 
n k~ P ~ n k = — P-V^n° = 5(fi — e k )p-v k , we obtain 

Jp = Yl 5nk ik = £ <*(M ~ £k)p-v k jk (2-15) 

k k 
= Yj 6 ^~ £ k)P- v k 



k 



Comparing eqs. ( |2.6| ) and (2.15), we obtain 

^-n(i) . (2.i7) 

dp v \m' J ^ 
On the other side, to estimate the change of the total energy 

5E P = £ k$n k + k ')S n k5n k *, (2.18) 

h h h> 
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the deviation 5n k has to be expanded up to the second order in p: 5n k = —J2i Pid ki n° + 
Y.i,jPiPjdk,dk j n° /2. Then we have 

+ f( k > k')5(fi - £k)S(fi - ev)p-v h v Vlx , 

(2.19) 

where the symmetry relation f(k, k') = f(k', k) was used in the second term. The first term is put 
into 

~ ( d k^k) p- V fc n° = J25(fi- e k )p-v k v k ^ (2.20) 

k k 

where we used J2k dku,{ £ kP'^ k n ) = 0- Hence we have 

Lp - X^G" - e k )p-v k v kfl 



k,k' 



As a result, from eqs. (2.16) and (2.21) we obtain 

dE„ 
dp 

and using eq. ( [2.17 ) we conclude 

d,J n „ 3 2 E, 



(2.21) 



•V = -XT', (2-22) 



n 



• (2-23) 



Qpv dp^dp 

As the assumption of an isotropic system, to regard f(k, k') just as a function of cos 6 = k • k'/kkf, 
is not used in the above derivation, the result fl2.23|) is generally applied even for an anisotropic 
system.0) 

§3. Fermi Liquid Theory II 

In order to describe the state \p) microscopically, we introduce the Green's function 

Gv(k) = r- (3-1) 

p u - e k + fj, + i sign(e fc _ p - n) 



We use the abbreviated notation such as k = (k,u) and k! = (fc',u/). In eq. (3.1), e k represents 
energy of the bare particle k. In the following, the energy of the quasiparticle k will be denoted by 
e* k . To calculate E p of the state \p), the Green's function (|3.1| ) is to be used instead of 

G°{k) 



to - e k + n + i sign(e fc - fj) 
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which is usually used to describe the ground state© In place of eq. (3.1), we may equivalently use 

G°Jk) = 1 — -. (3.2) 

p u- e k+p + n + 1 sign(e fc - /x) 

This Green's function describes the system which has the same Fermi surface as for the ground 
state, but with the shifted one-body energy £j~+ P instead of For example, the one-body part of 
the total energy of \p) is expressed in terms of G' p °(k) = Gp(k —p), 

d 4 k 



K = JX4*>! = - [ J j^ £ * g pW 



d 4 k 
£k+pG p (k) 



(2vr)4 
d 4 k 



e k + P G°(k) (3.3) 



(2vr) 4 k+p 

= ^£ k+ p{c\c k )°. 
k 

To evaluate dE p /dp, the p-dependence in the argument of G p (k — p) can be neglected by virtue 
of quasimomentum conservation at each interaction point, or we can replace G'p(k) in E p by Gp(k) 
without loss of generality. Physically, this is related to the fact that replacing the distribution 
function n k by n k ^ p in the expression for E p is equivalent to replacing e k by £k+ p - This point will 
be seen below in §|| for a simple case. As a result, for the energy correction AE p = E p — E® caused 
by electron-electron interaction, we obtain the equation 

dAE p _ . f d 4 k dG° p (k)^ 



where E' (k) is the total self energy including improper as well as proper parts. In terms of the 
identity 

dG°Jk) n/ x „ 

-J±-l = Vk+pilG ° p (k)\ (3.5) 

we have 



dAE„ r d A k 



p - -i 



dp P J (2vr) 

d 4 k 



v k+m G°Jk) 2 Z'(k) (3.6) 



= "i J j^v k+pfl G° p (k)X p (k)G p (k), (3.7) 

where we used the relation T,' p (k)Gp(k) = T lp (k)G p (k) for the proper self-energy X p (fc) and the 
dressed Green's function defined by 

Gpik)- 1 = G%k)~ 1 -^ p {k). (3.8) 

Thus, together with 

dE° r d 4 k 

P ~ ^ *rtW3S(fc), (3-9) 



dp„ J (2tt) 
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eq. (3.7) leads to 



dE p . f d 4 k 



dvy, J (2tt) 

where we used 



: v k+Pfl G p (k), (3.10) 



G p (k) = G°Jk) + G°Jk)Z p (k)Gp(k). (3.11) 



The total current of the state \p) is just given by eq. (3.10), since we find 

d 4 k 

,v k (c k c k ) p = -i / ■ 

k 



j p = ^ v k{c{c k ) P = -ij -p^jv k G' p {k) 



-^v k+p Gp{k). (3.12) 
Therefore we obtain 

< 3 - 13 > 

Now we derive the other results of the previous section on the basis of the microscopic Fermi 
liquid theory using well-known Fermi liquid relations. 

Eli) 

To calculate a linear response to the 

vector potential A, the total current operator is defined by 

J„ = v„-eJ2^u A u, (3-14) 

V 

where 

v^ = J2v kfl clc k , Vkfl = -^, (3.15) 



and 



Then the dynamical conductivity is given by 

e 2 (. 



= t4tK ( + iO) + h / tSi #SrG(fc) e ^ ] , (3.17) 



where K^ u (u), defined by 

i Z" 00 

K nAu) = ~n {0\Tvp(t)v v \0)exp(kjt)at, (3.18) 

corresponds to the uniform limit k — > of the current correlation function K^ v {k). For the Green's 
function G(k) of the ground state, we omit the subscript p (= 0) for simplicity. In terms of the 
vertex function A v (k'; k), K^ v {k) is given by 

K^{k) = -~ J ^jv klx G{k' + k/2)G(k' - k/2)K u {k'- k). (3.19) 
The vertex function satisfies the following equation; 

A lx (k';k)=v k ,„-i J ^-Y^{k'X)G{k" + k/2)G{k" - fc/2)A M (fc"; fc), (3.20) 
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where T^(k,k') is the irreducible four-point vertex function. As is well known, in the limit k = 
(k,iv) — > 0, which corresponds to the k-to region in the vicinity of the Fermi surface, the product 
of the Green's functions G(k' + k/2)G{k' — k/2) does not behave regularly but it is written 



G(k' + k/2)G(k' - k/2) = G{k'f + 2iri4> 



',2 , •■)_;.. 2 k ' V k' 

uj-k-vt, 



(3.21) 



The first term G(k') 2 represents a regular part of the left-hand side. In particular, we are interested 
in the w-limit (k/to = 0) and the fc-limit (k/to = oo) of eq. ( |3.21| ), for which we obtain the relation 



{G{k'fY - {G{k'f} k = 27ri4% - e* k ,)5{u;'). 



In eq. ( 3.21 ), v% represents the velocity of the quasiparticle k, 

Vk ~ dk ■ 

The quasiparticle energy e£ is given as a pole of the Green's function, 

G~ 1 {k,e% - n) = 4 - e fc - X(£;,4 - fi)=0, 



(3.22) 



(3.23) 



(3.24) 



where the many-body effect in the one-particle spectrum is embodied in the selfenergy T,(k,uj). 
The renormalization factor is defined by 

8E(k,io) 



1 



dto 



uj=0 



Differentiating eq. (3.24) with respect to k, we have 



vl 1 



dZ(k,to) 



duj 



uj=0 



Vk 



dZ(k,to) 



dk 



uj=0 



thus we find 

V*k = Zk [Vk + 

Next we consider the derivative of the selfenergy 

dE(k,u) 



8E(Jfe,0) 
dk 



dk L 



for which, using 



we obtain 



d 



dk 



7 G(k',u') = {G(k') 2 } k lv k 



i i dT,(k',uj') 



dk'^ 



8E(k,uj) 
dk,,. 



d 4 k' 
(2*)' 



: r(°)(jfe,fc'){G(*') 2 }* U 



dH(k',u') 

aiP 



(3.25) 

(3.26) 
(3.27) 

(3.28) 

(3.29) 
(3.30) 



Note that the fc-limit {G(k') 2 } appears in eq. fl3.29| ) when the derivative is taken with respect to 
k' for fixed to' . Comparing eq. ( |3.30D with 



A^(fe') = Vk'fi - i 



d 4 fc' 



-T^(k',k"){G(k"f} k A^k"), 



(3.31) 
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which is derived from eq. ( |3.20| ), we find 

, k ( ij 



dE(k',u' 
dk'„ 



hence, 



because of eq. Q3.27 ). 



(3.32) 



(3.33) 



Using the above results, we can write the fc-limit of K fll/ (k), eq. (3.19), as 



K, 





r d 4 k' 




1 (2tt) 4 ^ 




f d A k> 




1 (2^) 4 ^ 




d 4 k' dvk'fj, 


nj 


(2vr) 4 dk' v 



v k ^{G(k'f} k K k u {k') 

k Gik>) 
^G(k'), 



(3.34) 
(3.35) 
(3.36) 

where we used eqs. (|3.29| ) and (|3.32|) in the second line. Substituting eq. (3.36) for the second term 
in the parenthesis of eq. ( |3.17D , we can write the conductivity as 

o2 



ie 



' w + iO 

Thus, the effective mass m! in the Drude weight, 



(if^ + iO)-^ 



k 



Recx^w) = 7re ( — I 5{u) + 



0": 



mc 



is given by the difference between the w-limit and the A:-limit of K^ u (k); 



nr 



(3.37) 



(3.38) 



(3.39) 



In other words, the Drude weight is determined by the singular (coherent) part of G 2 (k'), eq. ( 3.22 ). 
To estimate eq. ( |3.39 ), from eq. ( 3.20| ) we obtain 

d A k' 



A£(fc) = A*(fc) - i / t^I 1 ^ fc ({g^) 2 }" - {G(fc') 2 } fe ) A*(tf). 



Thus, eqs. ( ^22|) , (|333| ), (EOOp lead to 



z k A^{k) = v * kfl + Y, f{k, k')v* klll 5(fi - 4') = .//■ 
fc' 



where 



Since we have 



f(k,k') = z k z k ,T"(k,k'). 



r uj =K k , _ 



d 4 k' 



nj (2vr) 4 



A*(*0 ({G(fc') 2 f -{G^) 2 }") A^(fc')) 



(3.40) 



(3.41) 



(3.42) 



(3.43) 
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which is derived from eq. ( |3.19| ), using eqs. ( 3.33; ), ( |3.22 ) and (3.41), for eq. ( p.39|) we finally obtain 

1 (3-44) 



in' 



in agreement with eq. (|2.6| ). Hence, following the argument given below eq. (|2.6| ) and using jt = 
Vk = k/m, we can conclude the absence of renormalization in the Drude weight in Galilean invariant 
systems. In particular, when dv^jdky = 5^ u /m, eq. (3.36) gives 



K 



1 f d k G(k') = ~8^. 



^ m nj (2vr) 4 
Therefore, in this case, the sum rule, 



ire 



Re<v(w)cL; = ^- (i"^ 



irne 



m 



J \M>1 



(3.45) 

(3.46) 
(3.47) 



is saturated with the Drude weight, hence K^ u = 0. In general, a lost weight in the coherent part 
is transfered to an incoherent part; 



n 



in' 



flu \"" / 

where we defined the effective mass in the total weight, 



1 f°° lmK u Ju') . 
— / f^-^-doj' > 0, 



00' 



[IV 



1 Jg, 

n \ ^ 



(3.48) 



(3.49) 



The incoherent part of the spectrum cannot be described in the framework of the Fermi liquid 
theory© 

In the same manner as we did for G(k), we can derive the Ward identity for the Green's function 
G p (k) defined in eq. (j3T§|). Let us introduce e* kp defined by 



£ kp ~ £ k+p — S p (fe, £ kp — fJ,) — 0. 

Differentiating this with respect to p, we find 



de 



kp 



dp 



p=0 



d_ 

dp~ 



z k\Vk+ — S p (fc,0) 



(3.50) 



(3.51) 



p=0, 



Similarly as in eqs. (HH) - (|3T33|) , using 



and 



dp. 



dpu P 



(k',u') = {G(k') 2 Y lv K 



+ 



dZ p (k',u;' 
dpu 



(3.52) 
(3.53) 
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we get 



d~E p (k',u;' 
dPu 



4:ijr 



d 4 £r 

(2tt) 



-T^(k',k"){G(k") 2 Y [v w , 



+ 



dZ p (k",uj") 



(3.54) 



Here the regular function {G(k") 2 } in the o;-limit should be used. This is because the pole e^ p of 
G p (k,uj) adiabatically moves to as a function of p without crossing the Fermi surface, since the 
Fermi surface for G p (k') is held fixed as a function of p. Thus, because of 

d A k" 



Vk' 



(2*)' 



-T^\k',k"){G(k") 2 \ K^{k h 



which is obtained from eq. ( |3,20 ), we find 



^(k')=v k> + 



and eq. ( |3.51| ) gives 



ds 



hp 



dp 



dE p (fc',q/) 
dpu 



Jkp,i 



(3.55) 



(3.56) 



(3.57) 



p=0 



where we used eqs. (3.41) and (3.51). Equation ( 3.57| ) is to be compared with eq. ( 3.23j ). In 
particular, under the assumption 



we have 



Vk + V k ' - V k+ g - Vk'-q = 0, 

dT, p (k,Lo) dT, p (k,Lo] 



+ 



0. 



(3.58) 



(3.59) 



du> dp 

It is easy to check this equation for a few diagrams of lower order. Therefore, in this case, we 
see from eqs. ( 3.51| ) and ( 3.57] ) that the quasiparticle current is not renormalized; = Vk- This 
equation led us to conclude the absence of renormalization in D. 
To derive eq. ( ggg ), we write the o;-limit of K^ u (k) as 

d 4 k' 



fjLV 



QJ (2tt) 4 
d 4 k' 



v h ^{G{k') 2 } A»(k') 
d 

Vk'a-^—Gp{k), 



QJ (2tt) 4 ^dp, 

which corresponds to eq. (3.35) for K™. Subtracting eq. (3.36) from eq. (3.60), 

d 4 fc' 



(3.60) 



[IV fJ,V 



nj (2^) 4 
i r d*k' d 



d_ 

^ dp v 



r^—GpW + ^Gik'] 



nj (2vr) 4 dp, 



■(vk'+p^G p (k') 



p=0 



Therefore, in terms of eqs. (3.12) and ( |3 . 1 3| ) , for eq. (3.39) we finally obtain 



n 



m 



1 dJ., 



[IV 



1 d 2 E„ 



O dp v Qdp^dp,, 



(3.61) 
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§4. Perturbation Theory 

4-1 Linear response theory 



The effective mass in the Drude weight, defined in eq. ( [2.6] ) or eq. ( 3.44 ), which is derived on the 
basis of the Fermi liquid theory, gives us a definite picture on how the Drude weight is renormalized 
by the many-body effect. However, to estimate the effective mass perturbatively, the formula (3.61) 
is practically useful. Not only to verify this formula but to see how Umklapp processes affect the 
Drude weight, we apply a finite-order perturbation theory in the following. 

In the spectral representation, the dynamical conductivity ( 3.17] ) is written as 

e 2 / i r°° 1 \ 

= -[n J_jO\T%(t)v v \0) expM)di - -(0|e^|0) . (4.1) 

For Vfj, and i" in this expression, we use eqs. ( |3.15 ) and ( 3.16j ), where the sum over the spin 
component a is implicitly assumed in the summation S^. Then, taking the limit uj — > 0, we get 

In this expression, the sum is taken over the exact excited states |r) of the system. The right-hand 
side of eq. (4.2) is derived also from the energy shift E p — Eq caused by the perturbation, 

Hp - H = ^2(e k+p - ek)c\ a c ka 

ku 

-p-^ + jE^C' ( 4 - 3 ) 



As a result we obtain the relation 



for 



E p = (p\H\p) = (0\H p \0), (4.5) 

where \p) is the ground state of H p . 

At this point, it is clear that the method given here has a close connection with Kohn's dis- 
cussion0i) that the effective mass defined in eq. ( |2.7D is derived by investigating variation of the 
ground state energy under changes in the boundary condition of a finite system. In a finite sys- 
tem the change of the boundary condition causes a uniform shift of discrete momentum quantum 
numbers. The point is that, in our theory, the chemical-potential shift eq. (4.3) should be treated 
as a perturbation. Since the perturbation adiabatically applied to the system cannot deform the 
Fermi surface, we cannot reach \p) from |0) by perturbation theory, but we are led to the excited 
state of the energy E p , eq. (|4.5Q , instead of the ground state of H p . In practice, if Eq (= (0|-ff|0)) is 
expressed in terms of bare quantities and interaction parameters, E p is easily obtained from Eq by 
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replacing the distribution function with nk- p - Thus, as a method to calculate the Drude weight, 
our theory is equivalent to Kohn's argument, so that eq. (3.61) is generally applied with the proviso 
that E p is regarded as the energy of the state \p) derived from |0). The investigation made in the 
preceding sections is still instructive, since it clarifies the physical meaning of the Drude weight D: 
We must assume D > for the Fermi liquid state to be stable against the shift of the Fermi surface 
in the momentum space. 

4-2 Drude weight 

To evaluate eq. ( |4.2| ) by perturbation theory, the ground state |0) is expanded up to terms of the 
second order in the interaction V' = V — (V) ^ 

|0) = |0) + |0)i + |0> 2) 

r^O h r-E 
|m V \r)(r\V'\ S )(s\V'\0) 
r^o( E r-Zo)(E s -E ) 

|0)o ^ \(r\V'\0)o\ 2 U7) 

2 ^ o (E r -E y 1 



To put it concretely, let us study the Hubbard model 

st s. 



h = t + v = J2 £ kcL d ka + J2 Uh ^i- ( 4 - 8 ) 

k,a i 



We assume that the effective mass in eq. (f4.2|) becomes isotropic in the x-y plane, m' = 
In the zeroth order in U, the second term of eq. (fl.2j ) gives 



^7)„=o(0|^|0)o (4.9) 





,2 



*5>w («») 



ka t 1 



(4.11) 



p=0 



where 



and 



or, similarly we have 



n dp 2 p 

n k = o(0|4 CT c fcCT |0) , (4.12) 

E P 0) = n k £ k+p = n k~p £ k, (4.13) 
kcr kcr 

n \ 1 ^-^ drik de^ 1 ^-^ dn^ 2 



m'J dk.dk, fi (414) 
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It is noted in eq. ( |4.13| ) that shifting the distribution n k — > n k _ p is formally equivalent to replacing 
Sk+p by e k for fixed n^. Accordingly, in eq. (4.10), (n/m')o is expressed as an average of d 2 e k /dk 2 
over the states below the Fermi level, while in eq. ( 4.14j) it is given by an average of U 2 . over the 
states at the Fermi level, as in eq. ( |3,44j ). To estimate l/m' numerically, eq. (4.10) is more suitable 
than eq. ( |4.14| ). 

For the correction {n/m')2 of order 0(U 2 ), from the first term of eq. (4.2) we obtain 



2 \(r\vx\0)i 



where 



En 



'fi3 



E 

fcl~fc4 



(1 - n kl )(l - n k2 )n k3 n k4 



(v klX + v k2X -v k3X -v kiX f8 kl+k2 _ k3 _ k4 , 

(4.15) 

(4.16) 



and from the second term of eq. 

^(i(0|Cl°)i + o(0|^|0) a + 2 (0|^|0) ) 

(1 - n kl )(l - n k2 )n k3 n ki ( d 2 e kl 



TJ 2 

— v 



a 2 e fc „ a 2 e fc , d 2 e k . 



As a result, 



and we find 



m' 



dk 2 



dk 2 



dk 2 



hi 



+ k 2-k3- k 4 



n 



rn 



eq. (P^D + eq. (4.17), 

i a 2 



.£(2) 

O dp 2 p 



p=0 



where 



U 2 



p n 2 J-i 

TJ 2 

n 2 z - 



(1 - n fcl )(l - n k2 )n k3 n ki 

^fcl+p ~T £ k 2+p ^/C3+p ^/C4+p 

(1 - nfe 1 _ p )(l - n k2 - p )n k3 _ p n k 



'kx+k%—kz—ki 



4-p 



fcl~fe4 



(4.17) 

(4.18) 
(4.19) 

(4.20) 
(4.21) 



The above results, following eqs. (4.11) and ( [4,19| ), verify the relation (3.61) up to terms of the 
second order in U. In particular, eqs. fl4.15| ) and (4.17) indicate (n/m')2 = for 



V k + V k > - V k+g - Vk'-q = 0, 



(4.22) 



i.e., the violation of the group- velocity conservation leads to renormalization of the Drude weight, 
in accordance with the note made below eq. ( 3.5§| ). Recently this point was discussed by Maebashi 
and FukuyamaEl) 

In Fig. 1, the ratio n/m' calculated with eqs. (4.11) and fl4.19|) are displayed as a function of 
n for U = and U = 1 for the square lattice with At = 1 and also for the case where the next- 
nearest-neighbor hopping t' = — 0.2i is included. In the figure, n/m for the total weight defined in 
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Fig. 1. (a) n/m' for U = 0, U = 1 and n/m for (7 = 1 are displayed as a function of n/2 for the Hubbard model 
in the square lattice with At — 1 (left), (b) The same quantities in the presence of the next-nearest-neighbor 
hopping t' — — 0.2i (right). The Drude weight D and the total weight of cr(u>) are proportional to n/m' and n/m, 
respectively. The curves for n/m' (U = 1) are not shown around the van Hove singularity. 



eq. (13191 ), which is given by eqs. (|4.14| ) and (4.17), is also shown. It is clear from the figure that 



ire 2 n 



the many-body effect modifies the Drude weight, 

D 

ml 

For U = 1, we see n/m' decreases precipitously as n approaches the van Hove singularity. In the 
second-order correction —{n/m') 2, we found that the part —(n/m')^ that is not due to Umklapp 
processes varies monotonically. For example, in Fig. 1 (a), — (n/m') 2 increases almost linearly up to 
less than 0.03 as n approaches half filling. The small and non-singular behavior of — (n/m!)\ due to 
normal processes is understood by observing that the relevant normal processes are not restricted 
in a narrow region around the Fermi level, where the effect of the van Hove singularity is prominent: 
For the Fermi surfaces of our example, energy conservation excludes normal processes at the Fermi 
level.® Thus in our example of a single-band model, the renormalization of the Drude weight 
and the conspicuous behavior shown in the figure are primarily caused by Umklapp processes, 
that is, the renormalization due to the normal processes violating the velocity conservation© are 
qualitatively unimportant here. 

We have dD/dn = ire 2 /m! > for the Drude weight in an isotropic system, since m'(= m) in 
this case is independent of n. Correspondingly, even in a lattice system, it is convenient to regard 
current-carrying carriers as particles (holes) when dD/dn is positive (negative). For instance, in 
the U = Hubbard model on a square lattice, we have dD/dn > for n < 1 as shown by the 
solid curve in Fig. 1, while dD/dn < for n > 1. Thus for n > 1 the total current may be well 
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regarded as being carried by holes in the filled band. On the other side, in one dimensional systems, 
D vanishes continuously not only for n — * and n — » 2 but also for n — * 1 for finite interaction 
U > 0.& Therefore, in this case, even in the underdoped regime n^l the current-carrying carrier 
is regarded as holes doped in an insulating state at n = 1. Prom this point of view, it is interesting 
to note that our result for n/m' (oc D) shows that dD/dn for U/At = 1 first decreases from a 
positive value as n increases from 0, and finally it becomes negative for n%,l. The result indicates 
a hole- like behavior around half filling, n^l. This behavior is caused by the fact that the slope 
dD/dn for [7 = vanishes as n — > 1, while a finite negative correction (n/m')2 decreases in this 
limit, as Umklapp processes becomes effective. When the Umklapp processes are effective, the 
Drude weight is reduced by the many-body effect, so that a large coupling constant p*U gives rise 
to a large difference (n/m')o — (n/m')2, where p* is the density of states at the Fermi level. Thus 
the decreasing derivative dD/dn for n^l as a function of U (> 0) is generally expected even in the 
weak coupling regime. For example, in our model, we assumed the coupling constant U/W = 0.5 
for the bandwidth W = 8t. It is noted however that our results do not imply D oc |1 — n\ for 
n — > 1 as in one dimensional systems, and that the hole-like behavior determined from the sign of 
dD/dn has nothing to do with the shape of the Fermi surface, namely, whether it is closed around 
k = (tt,tt). In effect, owing to the complete-nesting property occurring at n = 1, we should have 
to take account of antiferromagnetic ordering which sets in around half filling. Our results shown 
in Fig. 1 are valid when the ground state is a paramagnetic Fermi liquid. 

In a strong coupling regime, the decrease in the Drude weight around half filling is reproduced by 
a numerical technique© In our example of the weak coupling theory, however, the steep decrease 
in (n/m')2 of eq. ( |4.18| ) is mainly caused by eq. (4.15), as is clear from the figure. In other 
words, the behavior D — > as n — > 1 is ascribed to Umklapp processes becoming effective in this 
limit, rather than the decrease of the total weight n/m, eq. ( 3.4S| ). Therefore, in regard to the 
mechanism to reduce D, our result should not be simply compared with the situation in the strong 
coupling regime, as in the t-J model,© where the decreasing behavior of D would be mainly due 
to < D /ire 2 = n/m' < n/m — > 0. 

§5. Discussions 

From the above results, it is remarked that the large quasiparticle mass m* does not necessarily 
imply the enhanced mass m' in the Drude weight. This point is evident from the fact that, in an 
isotropic system, the former is affected by the many-body effect while the latter is not. Physically 
this means that a small Drude weight does not necessarily suggest a large specific heat coefficient. 
Thus, although m! and m* are related with each other (eqs. ( |2.13D and ( |2.14 )), they are practically 
regarded as independent quantities. In particular, this should be the case when we know little 
about the function f(k,k') by the first principle calculation. We introduced the mass m! for the 
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sake of convenience so that it represents the mass per particle of the system as a whole, 

E P P 2 



N 1ml 

This is obtained from eq. (3.61). In this form it may become even clear why electron-electron 
interactions do not modify ml in a Galilean invariant system; where the electron-electron interaction 
cannot change the inertial mass of the system as a whole. The effective mass ml thus must be 
positive for the system to be stable© As is clear from eq. (|2.13| ), to distinguish ml from m*, 
we must distinguish between the current and velocity of quasiparticle, and thereby we should not 
neglect the vertex correction. The difference is due to the Fermi liquid effect. 
If a Fermi liquid state were continuously driven to the metal-insulator transition, we should have 



ml — * oo. To see how this occurs, from eqs. ( |2.6j ) and (2.14) we may write 



' ' r i + %), (5-1) 



m' m* \ 3 
where 

^oE^'^d 1 -^). ( 5 - 2 ) 

m* 3 , 
k 

Ff = Ek,k' fe 'K • Vk>5(n - e k )5{n - e fc >) ^ 

3 J2k v k • v k 5(fi - e k ) 

Thus, to obtain ml — ► oo, we must have either (i) m* — > oo or (ii) F( — > —3. The former, if the 
coupling F( remains unaffected, corresponds to reduction of energy scale of the system, with no 
qualitative change of the low-energy physics. In this case, the total weight n/m, eq. ( |3.49| ), will 
be reduced as well. The latter, usually unnoticed, is a nontrivial possibility, which reminds us of 
the ferromagnetic instability Fft — ► — 1 in 3 He. In this case, the degeneracy temperature below 
which to validate the Fermi liquid theory will be severely suppressed as we approach the transition 
point D = 0. We must be in this situation to conclude D — > for finite n/m (> 0), and this 
was the case of our concern in the previous section. Nonetheless, in either case, the instability 
does not manifest itself as a thermodynamic phenomenon, for the vanishing of the Drude weight 
makes sense only in a coherent regime. Although the case (ii) and its physical consequences are 
interesting in its own right as a way to destabilize the Fermi liquid state, here we do not discuss 
this point any further but to point out the possible instability. In general, (i) and (ii) is to be 
regarded as independent possibilities, although these are related with each other in an isotropic 
system, for which m' = m < oo because of the Landau relation (1 + F(/3)/m* = 1/m. In a lattice 
system, on the other side, the relation does not hold so that in principle we can even think of the 
hypothetical system where the vertex correction due to Ff is neglected while the mass m* is heavily 
enhanced. In effect, the effective mass ml is generally anisotropic and the instability ml — > oo may 
be relevant particularly in an anisotropic system, e.g., in a quasi- two-dimensional system where 
Umklapp processes can become quite effective in the direction perpendicular to a two-dimensional 
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layer. Then, beyond the instability, if any, the Fermi surface has to be rearranged in the way the 
system avoids the instability. 

Finally we shall note that the plasma mode softens as uj p oc 1/ymf when m! — > oo: In a charged 
system, the dielectric function e{oS) is related to the dynamical conductivity, 

e(w) = 1 + —tr(u). (5.4) 



Therefore, eq. ( |2.5| ) gives 

Airne 2 



e M = - / 2 



from which the plasma frequency is given by 

4 = ^7-- ( 5 - 5 ) 



47me 2 



This result is valid only for u p <C £% , since we used the Fermi liquid formula for a(u). 

In summary, we derived a formula to calculate the Drude weight D on the basis of the Fermi 
liquid theory. To this end, we considered the state \p) which is obtained from the ground state 
by boosting it by p in the momentum space. Thereby we identified the instability of the Fermi 
liquid state caused by the metal-insulator 'transition' D — * 0. In a lattice system, electron-electron 
interactions enhance the effective mass m' , although an unrenormalized mass m! = m is concluded 
identically in a Galilean invariant system. As the enhancement is caused by Umklapp processes, 
the many-body effect on m' is largest in the vicinity of half filling n ~ 1 in a single-band model. 
This was shown for the Hubbard model in a square lattice by perturbation theory. As a result, even 
in a weak-coupling regime, we obtained the Drude weight D showing the behavior of the 'doped 
insulator' dD/dn < in the underdoped region n^l. It was noted that m' corresponds to the mass 
in the quasiparticle current jk and is generally different from the thermal mass m* defined through 
the velocity of quasiparticle v* k . 
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Appendix: Finite Temperature Formalism 

In this appendix, we outline the derivation of eq. ([[!]) by the finite temperature diagram tech- 
nique. The formalism used here is based on the work of Eliashberg. According to eq. (24) of 



18 



ref. 2C, the conductivity is given by 

f \ - ie2 Jv^ * 1 cosh- 2 (e* k - n/2T) „ 1^ ,1 cosh~ 2 (e|, — fi/2T)T 2 2(k, k'; uj) 
a ^ {u) ~ m j^^2T ^T2^* ^V+22.^2T ( w + 2i 7 *)(u; + 2i 7 *,) ^ 

(A-l) 

where 

T^-^ImEC^^Uo. (A-2) 
For our purpose we shall use only the real part of T22{k, k';u), for which eq. (12) of ref. gives 

ReT 22 (M» = (tanh £ * k ' ~ ^ + - - tanh ReT(k,k';u). (A-3) 



Therefore, noting 

and 7^ oc T 2 , in the collisionless regime 7^/0; —> of our concern, the limit u> — » leads to 



^cosh- 2 ^^^5(/,-4), (t-o) 



^(w) — - XI 1 - 4) + H /( fc > k ') v in v t'J{^ - 4M0 - 40 f ( A4 ) 

w i2 *; I it' J 

ie 2 1 

k 

where f(k, k') is defined in eq. ( |3.42| ). 

In the opposite hydrodynamic limit, w/7^ — > 0, the imaginary part of T22(k,k';u) plays an 
important role, and in contrast to the real part, to investigate ImT2 2 (k, k'; uj) poses a problem. 
This rather complicated task is indispensable, e.g., to evaluate the dc conductivity <r(0) (oc T -2 ) 
of a clean system at finite temperature© This point will be discussed in detail in the following 
paperS 
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